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Abstract 

We analyse the data for the proton structure function F2 over the entire domain, 
including especially low Q^, in terms of perturbative and non-perturbative QCD contribu- 
tions. The small distance configurations are given by perturbative QCD, while the large 
distance contributions are given by the vector dominance model and, for the higher mass 
qq states, by the additive quark approach. The interference between states of different 
qq mass (in the perturbative contribution) is found to play a crucial role in obtaining an 
excellent description of the data throughout the whole region, including photoproduc- 
tion. 



1. Introduction 

There now exist high precision deep inelastic ep scattering data 0, |^ covering both the low 
and high domains, as well as measurements of the photoproduction cross section. The 
interesting structure of these measurements, in particular the change in the behaviour of the 
cross section with at ~ 0.2GeV^, highlight the importance of obtaining a theoretical 
QCD description which smoothly links the non-perturbative and perturbative domains. 

In any QCD description of a 7*p collision, the first step is the conversion of the initial 
photon into a qq pair, which is then followed by the interaction of the pair with the target 
proton. Let a{s, Q"^) be the total cross section for the process 7*p X where is the 
virtuality of the photon and a/s is the 7*p centre-of-mass energy. It is related to th e forward 
7*p elastic amplitude A by the optical theorem, Im A = sa. We may write a double dispersion 
relation Q for A and obtain for fixed s 

-is, = E / J^^. I Jf^ Pis, M^ M-) I Im A,,,,is, M^ M^) (1) 

where M and M' are the invariant masses of the incoming and outgoing qq pair. The relation 
is shown schematically in Fig. 1. If we assume that forward qq + p scattering does not change 
the momentum of the quarksQ then Agg+p is proportional to (5(M^ — M'^), and (|I|) becomes 

/•oo <7/lf^ 

cris, Q') = E / 7W-o^2 Pis^ ag-g^,(s, M^) (2) 
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where the spectral function p(s, M^) is the density of qq states. 

Following Badelek and Kwiecinski 0] we may divide the integral into two partsQ, the region 
< Qq described by the vector meson dominance model (VDM) and the region > Qq 
described by perturbative QCD. Suppose that we assume pcgg^p is a constant independent of 
(which should be true modulo logarithmic QCD corrections) then the perturbative component 
of the integral is 

P- = / TTrrV-FvT—F^ P- = -is^Q +Qo)- (3) 



lQl{M^ + Q^y^ Jo (M2 + Q2 + g2)2 
Thus (0) becomes 

a{s, Q2) = a(VDM) + a^^D (s, + Ql) (4) 

where the QCD superscript indicates that the last contribution is to be calculated entirely from 
perturbative QCD. 



^In a more detailed treatment this assumption is no longer valid, see ( |2^ ) and (^) below, and the discussion 
in section 4. 

^Although Badelek and Kwiecinski base their fit to the data on eq. (||), they also discuss the more general 
case in which AI ^ AI' contributions may be included in the spectral function p. 
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We may use 
where x = Q'^/{s + Q'^ — M^) to rewrite (§) as 



^{s,Q') = -^F2{x,Q') (5) 



F2{x, Q') = F,{VDM) + — ^ {x, Q' + Ql) (6) 

where x = {Q"^ + Q^) /{s + Q"^ + Qq — M^). The vector meson dominance term has the formQ 

F,(VDM) - - X: (7) 

where My is the mass of vector meson V and where the sum is over the vector mesons which fall 
in the region <Ql. The vector meson-proton cross sections (Jv{s) can be determined from 
the np and Kp total cross sections using the additive quark model and jy from the leptonic 
width of the vector meson V. The last term in can be determined from perturbative 
QCD using the known parton distributions. This approach was first proposed by Badelek and 
Kwiecinski (BK) [Q. We see that the BK model, (§) and @, makes a parameter free prediction 
of F2{x,Q'^) which is expected to be valid, for s S> Q^, for all including very low The 
BK predictions give an excellent description of the F2 data for Q"^ > 1 GeV^, but overshoot the 
new measurements of F2 for smaller values of Q^. This deficiency of the model was removed 
in a fit to the F2 data performed by the HI collaboration but at the expense of using an 
unreasonably low value for Ql = 0.45GeV^ and of introducing an ad hoc factor of 0.77 to 
suppress the VDM term. 

The Badelek-Kwiecinski idea to separate perturbative and non-perturbative contributions is 
very attractive. To exploit it further we must achieve a better separation between the short and 
long distance contributions. To do this we take a two-dimensional integral over the longitudinal 
and transverse momentum components of the quark, rather than simply over the mass M of 
the qq pair. 

The contribution coming from the small mass region is pure VDM and is given by (^. 
However, the behaviour of the cross section at large is a more delicate question. The part 
which comes from large fc^ of the quark can be calculated by perturbative QCD in terms of the 
known parton distributions, whereas for small fcy we will use the additive quark model and the 
impulse approximation. That is only one quark interacts with the target and the quark-proton 
cross section is well approximated by one third of the proton-proton cross section. 

At this point it is interesting to note some recent excellent parametric fits of the data for 
F2, or rather for cr(7*p). The first is based on @ and the generalised VDM P]. To be more 
precise it is based entirely on a parametrization of the vector meson + proton cross section 
and does not take advantage of our present knowledge of perturbative QCD. As a consequence 



^Strictly speaking (Q) is the formula Ft- The small longitudinal component will be discussed later. 
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some anomalies appear. For instance the photoproduction cross section becomes negative for 
^/s < 6GeV (or cr(yp) < for My > 0.26-^5). Second the model has anomalously large values 
oi R = ai/aT (where Fl is obtained by including a factor ^Q"^ /My on the right-hand-side of 
formula (|^) for Ft)- In the well-known deep inelastic region the model predicts R > 1 for 
> 35GeV2 and x > 0.01 (and even R > 4 for x > 0.1) whereas the data indicate that 
R ~ 0.2 — 0.3. This effect probably reflects, as the authors note, the omission of allowing ^ to 
depend on Q^, see (|3^) below. Rather their model has ^ = 0.171 for all Q^. 

An earlier approach based on the generalised VDM can be found in ref. 0. In addition 
to the VDM contributions, this work contains a contribution at small x coming from "heavy" 
long-lived fluctuations of the incoming photon, which are parametrized in terms of a "hard" 
Pomeron whose intercept is found to be ap' = 1.289. 

Another fit of the Fa dat a is based on the Regge motivated ALLM parametrization 
The description, with 23 parameters, describes the data well and may be used to interpolate 
the measurements. On the other hand the physical basis of the parametrization is not clear. 
For example a variable xjp is defined by 

where W = ^/s is the 7*p centre-of-mass energy, M is the proton mass and Mjp reflects 
the energy scale of Pomeron exchange. This latter scale turns out to be extremely large, 
Mp = 49.5 GeV^, much larger than any hadron or glueball mass. Secondly the intercept, 
0^(0), of the secondary trajectory decreases with Q^, which is contrary to Regge theory (where 
aji is independent of Q^). 

The description of the F2 or o"(7*p) data presented in this paper is quite different. We use a 
physically motivated approach with very few free parameters, and we clearly separate the con- 
tributions to F2 coming from the large (small quark /ct) and small (large kx) distances. A recent 
study with a similar philosophy to ours can be found in ref. 0. They achieve a qualitative 
description of the experimental data over a wide range of photon virtualities (Q^) and energies 
{W) in terms of short and long distance contributions. They emphasize that even in the very 
low region the short distance contribution is not small, and also that at large the long 
distance effects still contribute. Here we present a quantitative study which involves a more 
precise approximation for the qq + p cross section and includes consideration of the longitudinal 
structure function Fl. Other differences are that we compute the (small kr) non-perturbative 
component using the VDM for small qq masses M < Qq and the additive quark model for 
M > Qo; we do not need an artificial suppression^ of the VDM component. Moreover we 
make a detailed fit to the F2 data in terms of an unintegrated gluon distribution which we de- 
termine using an unified evolution equation which embodies both DGLAP and BFKL evolution. 



*In ref. ||^ an ad hoc suppression factor of 0.6 is used. 
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2. The 7*j9 cross section 

The spectral function p occurring in (|l]) may be expressed in terms of the 7* qq matrix 
element Ai. We have p oc with, for transversely polarised photons, 



Jz{l - z) 
Mt = ^2 7^ ux{-fe±)uy 

(g± • fcr)[(l - 2z)\ T 1] ^A,-A' + Amg 5xx' 



(9) 



We use the notation of ref. IIQ], which was based on the earlier work of ref. [TT|. Namely the 
photon polarisation vectors are 

St = £± = ^ (0, 0, 1, ±0, (10) 

and A, A' = ± 1 corresponding to q,q helicities of ± |. Also we introduce 

= z{l-z)Q^ + ml (11) 

Note that (Q) is written in terms of "old-fashioned" time-ordered or light cone perturbation 
theory where both the q and q are on-mass- shell. This form is appropriate when discussing 
the dispersion relation (0) in the qq invariant mass. For high photon momentum the two 
time-ordered diagrams have a very different energy mismatch 

(aE 9l±Ml] « (Ai?' ~ p^) , (12) 



2p 
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and so the contribution from the diagram (A£") with the "wrong" time-ordering may be ne- 
glected. The remaining diagram, with energy denominator 1/AE, leads to the behaviour 
1/(Q^ + ^t) contained in (^, as can be seen on using ( [T4|) below. 

In terms of the quark momentum variables z, of Fig. 1, equations (|I|) and @ become 
= «TT E fd^ d'kr {MtM*t)N, - Im A^-^+p 

^ ~r7T \ _|_ -\ S 

(13) 



q ^" A = ±i 



where the number of colours A^^^ = 3, and Cg is the charge of the quark in units of e. We shall 
give the corresponding cross section cr^ for longitudinal polarised photons in the section 2.1. 
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The dispersion relation in has become, in (|T3[), a two dimensional integral. The 
relation between the variables is 

= (14) 

z{l-z) ^ ' 

where is the mass of the quark. For massless quarks z = \{1 + cos 6), where 6 is the angle 
of the outgoing quark with respect to the photon in the qq rest frame. The dz integration is 
implicit in (Q) as the integration over the quark angular distribution in the spectral function p. 

To determine F2{x, Q^) at low we have to evaluate the contributions to (Tt coming from 
the various kinematic domains. First the contribution from the perturbative domain with 
> Qq and large fcf^, and second from the non-perturbative or long-distance domains. 

2.1. The 7*p cross section in the perturbative domain 

We may begin with the two gluon exchange contribution to quark-quark scattering 

where ± It are the transverse momenta of the gluons. Thus for g-proton scattering we obtain 



(15) 



q+p 



J asm fix A) f- 



where 



/(x, /y) = a;9(7(a;, /|)/91n/^ 



(16) 



(17) 



is the unintegrated gluon density. The process is shown in Fig. 2. Finally for qq + proton 
scattering we have to include the graph for q+p scattering. For both the q and q interactions 
we have two diagrams of the type shown in Fig. 3 with A4*{kT + It) and A1(/ct)- We obtain 



E 



TT 



J d^kiT dz d^lT 



fix, I 



{l-zY + z' 



Iirp 



asilr) X 



' fciT It_ 



k 



IT 



Do 



1 



(18) 



where 



X = {Q^ + M^)/s, 



(19) 



Di = klj, + z{l - z)Q^ + 



9' 



L>2 = {It - kir f + z{l - z)Q^ + ml 



(20) 
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Expression ([T8|) is written as the square of the amphtude for quark- ant iquark production, 
where we integrate over the quark momentum ki^ in the inelastic intermediate state, see Fig. 2. 
The first term, proportional to I/-D1, corresponds to the amplitude where the gluon couples 
to the antiquark ^2, while in the second term, proportional to I/D2, the gluon couples to the 
quark ki. Of course form (^) can also be used to calculate the cross section for high kr dijet 
production (7*^ — > qqp), where kij- and fc2T refer to the transverse momenta of the outgoing 
quark jets. 

To separate the perturbative and non-perturbative contributions to the cross section (pTS] ) 
for our inclusive process we have to introduce a cut on the quark transverse momentum (as well 
as on the qq invariant mass M). At first sight it might appear that to obtain the perturbative 
contribution we simply require kiT > ko- However this implementation of the cut-off would 
not be correct. For instance if, as in Fig. 2, the two exchanged gluons couple to the ki line, 
then = It ~ ^it niay be small and in the limit — and small we would have an 
unphysical infrared singularity in the region of large kix and It-, but small /c2t, coming from the 



I/D2 term in (18). To see better the origin of the infrared singularities we perform the square 



and write the expression in curly brackets in ([T8|) in the form 



[(l-;^)2 + ^2]^^^.(^^_^^^)_^2 



(21) 



D1D2 



The danger comes from the second term, which corresponds to Fig. 2, whereas the last term, 
which describes interference, is infrared stable, as we will show later. Our aim is to separate 
off all the infrared contributions into the non-perturbative part. Therefore to evaluate the 
perturbative contribution coming from the second term we have to use the cut-off \lT—kiT\ > k^. 
This is equivalent to changing the variable of integration for the second term from kiT to 
It — kiT, and so its contribution is exactly equal to that of the first term. An alternative way 
to introduce the same cut-off is to separate off the incoming qq configurations with /c^ < /cq so 
that (|1^) becomes 



= V ^ / SkTdzdHT asit'r) 



(22) 



X 



[(1 _ + z^]kl + ml [(1 - zf + z^]kT ■ (kr + h) + mi 



(Q2 + A;2)2 (Q2 + ^2 )(g2 + (fc^ + f^) 



Note that the transverse momentum fc^ of the incoming quark is equal to fcij- when the gluon 
couples to the antiquark (first term in (pT])) and is equal to kix — It when the gluon couples 
to the quark (second term in (|2T|) ). Working in terms of the variable fcr corresponding to the 
dispersion cut shown in F ig. 1 has the advantage that it is then easy to introduce cut-offs with 
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respect to the invariant qq masses M and M', which we need to impose in order to separate off 
the non-perturbative VDM contribution^. 

Another argument in the favour of the cut written in terms of initial quark momenta kx 
comes from the impact parameter representation. Instead of kx we may use the transverse 
coordinate h and write the cross section (|22|) in the form 

ctt oc J dzd%\^^{b)\'^f{x,b)as{b) (23) 

where the gluon distribution 



/(:.,6)= /U^ll-e-'-"]/*!^. (24) 



The photon "wave function" is given by 

\^,ib)\' = Y.c^ey + il-zY]Q'K',{Qb), (25) 

where for simphcity we have set rriq = 0. The photon wave function is simply the Fourier 
transform of the matrix element A4 given by (^. It is most natural to take the infrared cut- 
off in coordinate space, say b < b^. The variable which is the Fourier conjugate of b is the 
incoming quark momentum knp of Fig. 1 (rather than the intermediate transverse momentum 
kiT of Fig. 2). This is further justification to impose the infrared cut in the form kf > fcg. 

Now let us consider the interference contribution, that is the last term in (^). It is infrared 
stable since in the limit and — * it takes the form 



/ 



d^krkT- (kT + h) _ f d{\kT + lT\)_ ^26) 



kUkr + ItY J h 



when \kT + It\ is small. We have used boundaries k'^ = k^ and = Qq to separate the 
perturbative QCD (pQCD), additive quark model (AQM) and vector meson dominance (VDM) 
contributions. As a result the 7*p cross section formulae, (P2|), is asymmetric between the 
ingoing and outgoing quarks. The origin of the asymmetry is the difference of the transverse 
momentum of the outgoing quark (fcy + l^) and the incoming quark (fc^) in Fig. 3. Such a 
graph therefore represents the interference between M and M' ^ M states. To obtain the pure 
pQCD contribution we require the incoming qq system to satisfy > Q\ and kx > k^. Ideally 
we would like to impose the same cuts on the outgoing qq system, namely 

M« = '"-^fl^r'"' > Ql (27) 



^Of course the use of the Feynman rules would yield the same result, but the time-ordered or light cone ap- 
proach with the incoming q and q on-shell is more convenient when we come to separate off the non-perturbative 



component in terms of fcr < and M, M' < Q, 



0- 
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and k'rp = {kx + It\ > ko- However in a small region of phase space, where It lies close to —k^, 
we may have M' < Qq and/or k'j, < k^. For this region we therefore have interference between 
the pQCD and VDM (or AQM) contributions. There is no double counting since neither our 
VDM or AQMQ components contain interference terms. This is fortunate because we cannot 
neglect the contribution from this small part of phase space of Fig. 3 without destroying gauge 
invariance, which is provided by the sum of the graphs in Figs. 2 and 3. We stress that the 
contribution coming from this limited region It close to — fc-p is infrared stable and hence it is 
small and has little impact on the overall fit to the data. 

So far we have only calculated 0"^. In the same way we may calculate the cross section for 
longitudinally polarised incident photons. In this case the relation analogous to (|1^) reads 



= 2^ ^ J dzdkT —2 — — (^q-q+p{kT), (28) 

which on evaluating dgg+p gives 

71 Jkl 1% 



(29) 



(g2 + klY (Q2 + fc2 )(g2 + (fc^ + 



iFiom. the formal point of view the integrals over and k"^ cover the interval to oo. For 
the It integration in the domain Zf. < /q ~ IGeV^ we may use the approximation 

as{ll)f{x,ll) = ^asill) f{x,ll). (30) 

For fcy < fcp we enter the long distance domain which we discuss next. To be precise we use the 
formula (^2]) and (|29|) to evaluate the cross sections only in the perturbative domain > Ql 
and fcy > /cq. We exclude the region < Qq and k^ > k^ from the perturbative domain as 
the point-like (short-distance) component of the vector meson wave function will be included 
in the VDM term. 

2.2. The 7*p cross section in the non-perturbative domain 

There are two different non-perturbative contributions. First for < Qq we use the 
conventional vector meson dominance formula (^ for Ft{x,Q'^). We also should include the 
longitudinal structure function Fl^x^Q"^). Fl is given by a formula just like (^ but with the 
introduction of an extra factor ^Q"^ /My on the right-hand side. ^(Q^) is a phenomenological 
function which should decrease with increasing Q^ . The data for p production indicate that 

^For the AQM contribution the interaction with the target proton is described by the forward elastic quark 
scattering amphtude and hence we have z' — k'j, = kr and M = M'. 
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^("^p) ^ 0.7 |]T3|, whereas at large the usual properties of deep inelastic scattering predict 
that 

^ . IS < M (31) 

So throughout the whole region the contribution of F^ is less than that of Ft- In order to 
calculate Fl (VDM) we insert the factor ^Q"^ /My in (|^) and use an interpolating formula for ^ 

with ,^0 = 0.7, which accommodates both the p meson results and the deep inelastic expectations 
of (PT]). However the recent p electroproduction, 7*p — > pp, measurements indicate that 
<^ l{p) / <yT{p) may tend to a constant value for large . We therefore also show the effect of 
calculating Fl (VDM) from (|^) using 

see Fig. 9 below. 

The second non-perturbative contribution covers the low /cy part of the > Ql domain, 
that is the region with < k^. Here we use the additive quark model and the impulse 
approximation to evaluate the aq-q+p cross sections in formulas ([T3|) and (pH]). 

2.3. Final formulae 

For completeness we list below the formulae that we use for the non-pQCD contributions 
coming from the /cy < /cq domain. When M < Qq, with — 1 — l.SGeV^, we use the vector 
meson dominance model. We have 

.,(VDM) = .a ^ E ^^Q.^^^ ?o (^J^ j (35) 
with ^0 = 0.7, see (0). For the vector meson-proton cross sections, we take 

O-p = O-a; = ^ [0"(7r+p) + Cr(7r"p)] 

= a{K+p)+a{K-p)-'^[a{n+p) + a{n-p)]. (36) 
Finally for M > Qo (and /cr < ko) we use the additive qua rk model and impulse approximation 

a.(AQM) = « E I / dzdk^T ^ ^'.fp^^.'' ^-^^ (^') (^7) 
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a,(AQM) = c.Y.§^ I dzdkl ^^^f_[^ f N^^,-,^r> {W') (38) 



(g2 + 4 

where for (Jqq+p we take, for the hght quarks, 

2 



a,-,+p{W') = -a,T,{s = lW^). (39) 

The "photon" wave function contains propagators hke l/((5 + k'^) and in impact parameter 
hx space it receives contributions from the whole of the hx plane extending out to infinity. 
On the other hand confinement restricts the quarks to have limited separation, say hx = 
\biT — 62t| ^ 1 fm. To allow for this effect we have replaced Q hj = Q + jj?' in (p7|) and 



|), where /i is typically the inverse pion radius. We therefore take jj? = O.lGeV^. This change 
has no effect for Q"^ ^ /i^ but for Q'^ < it gives some suppression of the AQM contribution. 

2.4. The quark mass 

In the perturbative QCD domain we use the (small) current quark mass mcurr? while for 
the long distance contributions it is more natural to use the constituent quark mass Mq. To 
provide a smooth transition between these values (in both the AQM and perturbative QCD 
domains) we take the running mass obtained from a QCD-motivated model of the spontaneous 
chiral symmetry breaking in the instanton vacuum [l^ 



The parameter A = = 1.09 GeV, where p = 1/(0.6 GeV) is the typical size of 

the instanton. fi is the natural scale of the problem, that is /i^ = z{l — z)Q'^ + or 
/i^ = z{l — z)Q'^ + [It + as appropriate. For constituent and current quark masses we take 
Mq = 0.35 GeV and m^urr = for the u and d quarks, and Mq = 0.5 GeV and mcu„ = 0.15 GeV 
for the s quarks. 



3. The description of the data for F2 

Though in principle it would appear that we have a parameter-freeQ prediction of -^2(3^, Q^) 
at low Q^, in practise we have to fix the values of the parameters fcg and Qq. Recall that k^ = k^ 
specifies the boundary between the non-perturbative and perturbative QCD components, and 
that = Qq specifies the boundary between the VDM and AQM contributions to the non- 
perturbative component. The results that we present correspond to the choice Ql = 1.5 GeV^, 
for which the VDM contribution is computed from the p, u and meson contributions (with 
mass My < Qq). The more sensitive parameter is /cq. We therefore present results for two 
choices, namely /cq = 0.2 and 0.5 GeV^, which show some interesting and observable differences. 
The results are much more stable to the increase of /cq from 0.5 to 1 GeV^. 

^Apart of course from the form of the input gluon distribution, g{x, Iq). 
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To calculate the perturbative contributions we need to know the unintegrated gluon dis- 
tribution /(x, /|.), see ( p2D and (pQ]). To determine we carry out the full programme 
described in detail in Ref. We solve a "unified" equation for f{x,l^) which incorporates^ 
BFKL and DGLAP evolution on an equal footing, and al lows the description of both small 
and large x data. To be precise we solve a coupled pair of integral equations for the gluon and 
sea quark distributions, as well as allowing for the effects of valence quarks. As in Ref. [0 we 
take 1^ = 1 GeV^, but due to the large anomalous dimension of the gluon the results are quite 
insensitive to the choice of Iq in the interval 0.8-1.5 GeV. 

The starting distributions for the evolution are specified in terms of three parameters N, A 
and j3 of the gluon 

xg{x,ll) = Nx-^{l~xf (41) 

where 1^ = 1 GeV^. At small x the gluon drives the sea quark distribution. The /ct factorization 
theorem gives 

S,{x, Q') = £ y / ^ SLA^, k\ Q^) f (|, Q^) (42) 

where S'box describes the quark box (and crossed box) contribution. The full expression for S'box 
is given in Ref. |]16| . Thus the sea Sq is given in terms of the gluon / except for the contribution 
from the non-perturbative region fc^ < k^, where we take 

S^P = 57 = 2S^^ = Cx-°-°^(l-x)^ (43) 

The parameter C is fixed by the momentum sum rule in terms of the parameters A^, A and f3 
specifying the gluon. The charm component of the sea is obtained entirely from perturbative 
QCD (see |jl6|) with the charm mass rric = 1.4 GeV. The valence quark contribution plays a 
very minor role in our analysis and so we take it from the GRV setQ of partons [IS]. Of course 



the sea quark distributions Sqi^x^Q"^) of (^) (and (0)) are used only to get a more precise 
determination of /(x, /^) through the coupled evolution equations. These forms for Sq are not 
used in our fit to the F2 data since the sea contribution is already embedded in (^) and (pQ]) . 

We determine the parameters A^, A and (3 by fitting to the available data for F2 with x < 0.05. 
We present two fits corresponding to a larger perturbative QCD contribution (Fit A with 
kl = 0.2 GeV^) and a smaller pQCD component (Fit B with /cq = 0.5 GeV^). The values of 
the gluon parameters are given in Table I and the quality of the description of the F2 data is 
shown in Fig. 4. Only a selection of the data fitted are shown in Fig. 4. Both descriptions 
are in general satisfactory, but Fit A is superior mainly due to Fit B lying below the data for 

~ 1 GeV^. This difference is better seen in Fig. 5 which shows the fit as a function of for 
various fixed values of x. We see that Fit A, with the larger perturbative component is more 



^Following Ref. pq | we appropriately constrain the transverse momenta of the emitted gluons along the 
BFKL ladder. There is an indication, from comparing the size of the next-to-leading ln(l/a;) contribution 
to the BFKL intercept with the effect due to the kinematic constraint that the incorporation of the 
constraint into the evolution analysis gives a major part of the subleading ln(l/a;) corrections. 

^The GRV valence distributions were fitted to the MRS(A) distributions ^ a.t = 4 GeV^. 



11 



Table 1: The values of the gluon parameters of eq. (^T]). 



(GeV2) 


N \ P x^/datapoint 
[423 points] 


Fit A 0.2 
Fit B 0.5 


0.97 0.16 3.6 1.09 
0.42 0.32 3.7 1.70 



able to accommodate the charge in slope going from high to low Q^. It is informative to show 
the components of the cross section. The breakdown is shown in Figs. 6 and 7 for Fits A and 
B respectively for the maximum energy W = 245 GeV for which data are available. It appears 
that the low behaviour of the pQCD component with low 1^ plays the vital role. 

The description of the F2 data by Fit A is better than that obtained by Badelek and 
Kwiecinski which is to be expected since we perform a fit to the data, albeit with a very 
economical parametrization. Fig. 5 also shows the HERA photoproduction measurements at 
W = 170 and 210 GeV. These data are not included in the fit. We see that our description 
overshoots the published HI ||2l| and ZEUS p2| measurements, although by a smaller margin 



than that of ref. [Q. On the other hand our extrapolation is in excellent agreement with a 
subsequent analysis of ZEUS data performed in ref. |2^. We will return to the comparison 
with photoproduction data when we study the effects of a different choice of the quark mass. 

4. Discussion 

We have made what appears to be in principle a prediction of F2, or rather of a^*p, over 



'■y'p, 

12 



the entire Q range which relies only on the form of the initial gluon distribution, see (^) and 
the parameter values of Table 1. However a comparison of the results of Fits A and B show 
that in practice the results are dependent on the choice of the boundary = between the 
perturbative and non-perturbative contributions, where ifc-r are the transverse momenta of 
the incoming q and q which result from the 7* qq transition. 

There are compelling reasons to select Fit A with k^ = 0.2 GeV^, which has the larger 
perturbative QCD contribution. Fit A is not only preferred by the data, but it also yields an 
input gluon with a more reasonable small x behaviour. In fact for Fit A {k^ = 0.2 GeV^) the 
AQM contribution is almost negligible and the fit produces a reasonable A, namely A = 0.16. On 
the other hand Fit B (with k^ = 0.5 GeV^) requires a larger A, A = 0.32, in order to compensate 
for the much more flat behaviour of the rather large AQM component. Further support 

for Fit A comes from the predictions for the longitudinal structure function, Fl. Fig. 8 shows 
that the prediction from Fit B is much larger than that of Fit A due mainly to the large AQM 



contribution. Fig. 8 also shows the expectations for F^ from the analysis of ref. p4| and from 



the MRST partons |25] of the most recent global parton analysis. We see these independent 



determinations of Fl favour the prediction of Fit A. 
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For completeness we show by the dashed curve in Fig. 9 the predictions of ai/ax versus 
obtained from Fit A. This figure also shows the effect of replacing (|32| ) by (p3| ) in the formula 
for the VDM contribution to F^. Recall that ( p3| ) was motivated by the possibility that the 
ratio <7l{p)/o'l{p) for p meson electroproduction tends to a constant value A as Q'^ ^ oo. We 
see from Fig. 9 that this change to the VDM contribution affects F^, and hence oi^joT^ mainly 
in the interval 0.2 < < 10 GeV^. It is straightforward to deduce from Fig. 9 the effect of 
changing the value of the parameter of (^3|) to match the constant limit A observed for the 
p ratio. 

A remarkable feature of the recent measurements [|l], |^ of o"(7*p) = (47r^a/(5^) F2{x,Q'^) 
at fixed W, is the transition from a flat behaviour in the low domain to the steep ~ 
fall off characteristic of perturbative QCD, see Fig. 5. The transition appears to occur at 

~ 0.2GeV^. Such a break with decreasing may reflect either the saturation due to the 

onset of absorption corrections or the fact that we are entering the confinement domain. The 

observed features of the data favour the last possibility. First there is no similar break in the 

behaviour of F2 as a function of x at low x which would be expected if absorptive corrections 

were important. A related observation is that the break, as a function of Q^, appears to occur 

at the same value ~ 0.2GeV^ for those W values for which data are available. Moreover we 

directly estimated the effect of the absorptive corrections using the eikonal rescattering model 

and found that they give a negligibly small effect on the behaviour of the cross section and 

of F2. On the other hand, if the break is due to confinement then it is expected to occur at a 
2 

value of Q which corresponds to the distances of the order of 1 fm, that is 

z{l-z)Q^ ~ gV5 ~ (0.2GeV)2 (44) 

which gives ~ 0.2GeV^ where the break is observed. 

In our calculations we have used a running quark mass which links the current (mcurr) to 
the constituent (Mq) mass. The growth of rrig in the transition region from perturbative QCD 
to the large distance domain is an important non-perturbative effect, which we find is required 
by the F2 data. /^From the theoretical point of view such a behaviour of rrig may be generated 



by the spontaneous breakdown of chiral symmetry in the instanton QCD vacuum ||T5[. The 



qualitative features are that rriq ~ Mq if the virtuality q of the quark is less or of the order of 
the square of the inverse of the instanton size, but that rriq decreases quickly as increases. 
In our analysis we have used a simplified power approximation for m^, see (^Oj). 

It is interesting to explore the effect of a different choice of quark mass. The dashed curves 
in Fig. 10 show the effect of using the constituent (fixed) mass Mq of the quarks in all the 
contributions to F2 or (j(7*p). As expected in the large 3> Mq perturbative domain the 
change has little effect. For small it reduces the predictions. For example, the photop roduc- 
tion estimates for W ~ 200 GeV are reduced by more than 10% and would bring our analysis 
more into line with the published HI and ZEUS photoproduction measurements. However our 
running quark mass predictions (continuous curves) are more physically motivated and should 
be more reliable. It will be interesting to see if their agreement with the experimental values 
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extracted in ref. is maintained when the new photoproduction measurements are available 
from the HERA experiments. 

A noteworthy point of our description of the F2 data is the importance of the non-diagonal 
(M 7^ M') perturbative QCD contribution to the double dispersion relation (|l]). The contribu- 
tion, which comes from the interference terms in (|22|) (and (p9D ), corresponds to the diagram 
shown in Fig. 3. It clearly has a negative sign, and moreover 

After the integration over the azimuthal angle in (P^D, the interference term exactly cancels 
the diagonal first term for any < in the limit of — and nig = 0. As a result 
the perturbative component of the cross section coming from the region of small It essentially 
vanishes[3 as 0. This property, seen in the It < lo components shown in Figs. 6 and 7, 

helps to reproduce the very flat behaviour of cr(7*p) observed at low Q^, Q"^ < 0.2GeV^. The 
fact that this low It gluon contribution becomes very small as decreases (and in fact vanishes 
for It < kT in the ^ limit) may be considered as a justification of the perturbative QCD 
contribution to F2 for low Q^. The VDM cross section (and other diagonal contributions as 
well) decrease as l/(My -f- Q'^Y require just such a component which increases with in 

order to compensate the decrease of the diagonal terms. The compensation is well illustrated 
by Figs. 6 and 7 which show the behaviour of the various components as a function of Q^. 
Of course the compensation (that is the effect of the vanishing of the low contribution as 
0) is more manifest in the Fit A where a larger part of the phase space is described in 
terms of perturbative QCD. 

It is interesting to note that in this paper we have included two different types of interference 
effect. First we have the dominant interference between the large M and M' states which 
gives rise to the decrease of the pure perturbative small It component of the cross section as 

— * 0, and which is responsible for the good description of the low data. Then there 
is the interference between the perturbative and non-perturbative amplitudes which we have 
modelled using the perturbative formula in the region of small M' and/or small \kT + It\- We 
have noted that this contribution is small due to the infrared stability of the integral, as was 
shown in (p6|). 

In summary we obtain an excellent description of F2, or rather of a^*p, over the entire 
range (from very low to high values of Q^) in terms of physically motivated perturbative and 
non-perturbative contributions. The choice of the boundary between the perturbative and non- 
perturbative domains which gives an excellent fit to the data, is also found to yield a sensible 
gluon distribution and reasonable predictions for Ft- 



^^Oi course there is also a non-negligible contribution coming from the domain It > kT which does not vanish 
as g2 ^ 0. 
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Figure Captions 



Fig. 1 The schematic representation of the double dispersion (|l|) for the 7*p total cross section 
a{s, Q^) at fixed cm. energy ^/s. The cut variables, M and M', are the invariant masses 
of the incoming and outgoing qq states in the quasi-elastic forward amplitude, Agg^p. 

Fig. 2 The quark-proton interaction via two gluon exchange. The spectator (anti)quark is shown 
by the dashed line, /(x, l"^) is the unintegrated gluon distribution of the proton. 

Fig. 3 A "non-diagonal" qq — proton interaction. 

Fig. 4 The description of the F2 data obtained in Fits A and B. Only a subset of the data fitted 
is shown. 

Fig. 5 The curves are the values of the virtual photon-proton cross section cr^.p of as a 
function of for various values of the energy W = ^/s corresponding to Fits A and B 
(multiplied by the factor shown in brackets). The data |^ are assigned to the value of 
W which is closest to the experimental W bin. The upper, lower photoproduction (solid 
triangular) data points correspond to W = 210, 170 GeV and are from the HI and 



ZEUS [E2| collaborations respectively. The open triangular points are obtained from an 



analysis of ZEUS photoproduction data reported in a thesis by Mainusch p3|. 



Fig. 6 The various components of defined in Section 2.3) shown as a function of at 

W = 245 GeV for Fit A (with /eg = 0.2 GeV^). The bold curve shows their sum, a^.p, 
compared to the HERA measurements. 

Fig. 7 The same as Fig. 6 but for Fit B (with k"^ = 0.5 GeV^). Th e poorer description of the 
data in the region ~ 1 GeV^, as compared to Fit A, is clearly apparent and can be 
attributed to the smaller perturbative QCD component at low gluon l^. 

Fig. 8 The predictions for versus ai W = 210 GeV from Fits A, B (with fcg = 0.2 and 
0.5 GeV^ respectively), together with the values obtained by Badelek, Kwiecinski and 
Stasto PH and from the MRST set of partons p5|. 



Fig. 9 The dashed curve is the prediction for a^/ax versus Q"^ at W = 210 GeV from Fit A. For 
comparison the continuous curve is the prediction obtained using a different choice of the 



VDM contribution to Fl; namely using (33) in the place of (^ 



Fig. 10 The dotted curves show the effect of using a (fixed) constituent mass, Mq, in all contri- 
butions. The running mass fit (continuous curves) and the data are those of Fig. 5. 
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